Determination of the (3x3)-Sn/Ge(lll) structure by photoelectron diffraction. 
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I. INTRODUCTION 

At a coverage of about 1/3 monolayer (ML), Sn and Pb 
deposited on Ge(lll) below ~ 550 K form a metastable 
(VS X a/3)R.30° phase, where the Sn or Pb atoms occupy 
T4 sites above the Ge lattice. This phase continuously 
and reversibly transforms into a (3 x 3) one, upon cool- 
ing below ^ 200 K. This new phase was first observed 
by scanning tunneling microscopy (STM) on the Pb/Ge 
and Sn/Ge systems and attributed to the. manifestation 
of a surface charge density wave (CDW) The nature 
of this phase transition is still an open issue, also because 
of the controversaLdetermination of the atomic structure 
of the two phases.Qcl The low temperature (3 x 3) phase 
is formed by three inequivalent Sn atoms per unit cell. 
X-ray diffraction (XRD) and low energy electron diffrac- 
tion (LEED) experiments now consistently point to an 
atomic structure where rOae Sn atom out of three pro- 
trudes above the surface.lai3 A vertical ripple of ~ 0.3 A is 
thus observed between the Sn atoms of the surface unit 
cell. Within the CDW model the inequivalence between 
the Sn atoms should disappear in the room temperature 
(RT) (VS X \/3) phase. This is consistent with the most 
recent experiments that indicate equivalent Sn ajtnms at 
the same height level in the adatom layer at RTJj'Q 

X-ray photoemission spectroscopy (XPS) experiments 
on the (3 x 3) phase indeed show two Sn ^ci core level 
doublets split by ~ 0.4 eV, with an intensity ratio of 
approximately 2:1 between the two components (A rand. 
B, majority and minority components, respectively) , El t2l 
thus indicating two different kinds of Sn atoms. This 
spectrum remainsjppactically unchanged throughout the 
phase transition.Ercl'U This opens the way tp-|alternative 
explanations, such as dynamical fluctuationsllil and order- 
disorder phase transitions ,□ where the local structure of 
the two phases is the same. Very recently the strong 
influence on the phase transition-oLsmall defect concen- 
tration has been put in evidence Ji3'c3 Recently published 
XPS studiesEj suggest that, at exactly 1/3 ML coverage, 
the Sn spectrum is different from that reported in the 
literature up to now. In particular the spectrum would 
show an intensity ratio between the A and B compo- 
nent appreciably different from that reported in previous 
works and would have an additional component in the 
high temperature {^/S x \/3) phase. These findings were 
used to support the CDW model. 

It is therefore important to determine the exact line- 
shape of the Sn Ad spectrum and to measure the (3 x 3) 
structure by XPS experiments on the single Sn core level 
components in order to check its consistency with that 
obtained by XRD and to identify the bonding sites of the 
atoms contributing to the majority and minority compo- 
nents. We have thus performed photoelectron diffrac- 
tion (PED) measurements of the Sn core levels at the 
ALOIS A beamline (Elettra Synchrotron). The wide va- 
riety of scattering geometries accessible at the ALOISA 
end station allowed us to measure separately the bond di- 



rection and length between the Sn adatom and its three 
nearest neighbour Ge atoms (Ge„„) and the Sn vertical 
distance to the next nearest neighbour Ge atom beneath 
(Ge„„„). In a previous paper we qualitatively reported 
on the preliminary experimental results, i.e. the substan- 
tial equivalence of the Sn-Ge„„ bond angle and length 
for both A and B type Sn atoma^ which mainly differ 
for their vertical height Sn-Ge„„„.E3 Here we present the 
quantitative determination of the bonding geometry of 
Sn atoms in the (3 x 3) unit cell, as obtained, by fitting 
to multiple scattering calculations (MSCD)^ the PED 
data obtained in a more extended energy range and with 
a better statistics. We attribute each of the two compo- 
nents of the Sn 4d photoemission spectrum to a specific 
bonding geometry of the two inequivalent types of Sn 
atoms. To our knowledge, such a direct correlation was 
never extracted before. The rippled structure obtained 
in the present work is in SCjupd agreement both with previ- 
ous XRD determinationsQu and with recent density func- 
tional calculations (DFT) performed in the local cLensity 
approximation with gradient correction (LDA-GC)EJ asd 
in combination with local-orbital (LDA-LO) methodslia. 

II. EXPERIMENTAL SET-UP AND SAMPLE 
PREPARATION 

The present experiment has been performed, in the 
UHV end station of the ALOISA beamHnell3'e3 at the 
Elettra synchrotron (Trieste). Details about the expeiii 
mental setup and sample cleaning are given elsewhere.E£l 
We have used a Reflection High Energy Electron Diffrac- 
tion (RHEED) system to monitor in real time the prepa- 
ration of the a-phase of Sn on Ge(lll). We have followed 
a procedure slightly different from the standard one re- 
ported in the literature. Sn is usually evaporated on sam- 
ples held at RT and successively annealed to T^ ~ 500 K. 
In this conditions, only a diffuse diffraction pattern is 
observed during deposition and the surface quality can 
be only judged a posteriori. Therefore the same surface 
coverage may not be easily reproduced, being subject to 
any unstability of the evaporation cell. By He atom scat- 
tering (HAS) experiments with high angular resolution, 
we have seen that the best surface quality (i.e. narrow- 
est diffraction peaks of both the (3 x 3) and x a/3) 
phases) is obtained when depositing on the surface kept 
close to 500 K.B At the ALOISA end-station we have 
used the RHEED system (which is a long range order 
probe similar to HAS) to monitor the diffraction pattern 
along the < 112 > surface direction while depositing at 
Ts ~ 500 K. The c(2 x 8) spots-|Soon disappear leaving 
only a streaky (2 x 2) pattern.E3 After 2/3 of the total 
exposure used in our work, the spots of the {\/3 x \/3) 
start to appear (see panel b of Fig. |^), but the deposition 
is just stopped when the half-integer streaky peaks disap- 
pear (see panel c of Fig. |^). The surface is shortly flashed 
at about 530 K and then cooled down. Higher flashing 
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temperature irreversibly leads to a stable (7 x 7) phase, 
as previously reportedcJ. Shorter exposures never lead to 
the complete disappearance of the (2 x 2) pattern. After 
this procedure sharp spots of the (\/3 x -\/3)R.30° phase 
are observed along the < 112 > direction. To detect the 
spots characteristic of the (3x3) phase, the RHEED pat- 
tern was observed along the < 110 > direction, as shown 
in the d panel of Fig. |l|, which is taken at T^. ~ 150 K. 
We have seen that no degradation of the photoemission 
spectra nor of the RHEED patterns occurs for at least 
10 hours, when the sample is left inside the experimental 
chamber. In any case, a new surface was always prepared 
after 10 hours. 

The whole set of spectroscopical data hereafter shown 
has been taken with constant electron energy resolution 
of 170 meV. The photon energy resolution has been al- 
ways < 150 meV. The overall instrumental resolution is 
thus < 225 meV. 



III. RESULTS AND DISCUSSION 

Fig. H shows a typical Sn 4d photoemission spectrum 
taken at ~ 210 eV photon energy on the, Q,x 3) phase. 
In agreement with previous observations,Bilil but.in con- 
trast with the results presented by Kidd et al.jlil it can 
be fitted to two spin-orbit split doublets, shifted by 
0.38 ± 0.02 eV and with an intensity ratio close to 2:1. 
Any deviation from this shape was always associated to 
residual (2 x 2) or (7 x 7) phases as observed by RHEED 
patterns. A spectrum taken at a deposition stage close 
to that shown in the panel b of Fig. can also be seen in 
Fig. H This spectrum shows a 4^5/2 peak at 24.3 eV 
binding energy, which is characteristic of the (2 x 2) 
phase. This observation is fully consistent witli-tlie. pre- 
vious measurements reported in the literaturclallO'til The 
experiments have been performed by taking the whole 
Sn 4d spectrum for each point of the FED angular and 
photon energy scans at Ts ^ 150 K. Four gaussians have 
been used to fit the spectra, while both the spin-orbit 
splitting of the doublets and the energy splitting between 
the A and B components were kept fixed. The same 
full width at half maximum has been used for all the 
components. A low order polinomial (n < 3) has been 
used for background subtraction. The FED data points 
shown in-the present analysis represent the so called x~ 
function jHJ x — , where / is the integrated intensity 
of the 4^5/2 component of the A or B type Sn atoms. Iq 
is a smoothly varying background of /, modulating both 
the angular and energy FED features. Iq essentially takes 
into account the energy dependence of the atomic pho- 
toionization cross-section and the angular dependence of 
the polarization and of the illuminated sample area. Iq 
has been evaluated by a low order (< 5) polynomial fit- 
ting of /. 

First we have measured the polar scans from the sur- 
face normal (normal emission, NE, conditions) to the 



horizon (grazing emission) for six different photon ener- 
gies and the three inequivalent main symmetry directions 
of the surface. The scans have been performed by keeping 
the sample in transverse-magnetic (TM) polarization, i.e. 
with the scattering plane normal to the magnetic field of 
the photon beam, and rotating the analyzer in the scat- 
tering plane. When the surface is oriented with the |TT2] 
direction in the scattering plane, at a specific angle a 
from the surface normal the analyzer is aligned with the 
Sn-Ge„„ bond angle (see Fig. |3|). In this, so called, bond 
emission (BE) geometry, the Sn photoelectrons can be ef- 
ficiently backscattered by the Ge„„ atom and give rise to 
intensity maxima and minima, when changing the photon 
energy. By a simple visual inspection, the bond direction 
a for both A and B Sn atoms can be roughly estimated as 
50° from the surface normal (see the polar scans taken at 
360 and 220 eV in Fig. ^). Given the information on the 
bond direction, the bond length I has been determined 
by variable energy FED scans in BE geometry. The po- 
larization vector of the photon beam has been oriented 
at 50° from the surface normal, along the bond direction. 
This choice enhances the sensitivity to the-Sn-Ge„„ bond 
length, as shown by MSCD simulations.t^ A very simi- 
lar photon energy dependence has been found for both A 
and B components in BE condition (see upper panel of 
Fig. H) .Ej Finally, the Sn height d above the underneath 
Ge„„„ atom has been measured by taking energy scans 
in NE conditions and with the sample surface in TM po- 
larization. In this case the A and B components display 
a remarkably different energy dependence, as can be seen 
in the lower panel of Fig. ^. 

The quantitative evaluation of I, d and a for both 
A and B type of Sn atoms has been performed by fit- 
ting the data to the MSCD simulations, starting from 
the geometric, structure recently obtained by LDA-GC 
calculations .E-l The simulations have been performed on 
a 30 atom cluster centered on the Sn atom by consider- 
ing multiple scattering up to the 8th order, 4tli Rehr- 
Albers order and with a pathcut of lxl0^^£3 Lower 
pathcut values did not modify the simulations apprecia- 
bly. A Debye temperature of 230 K and an inner poten- 
tial Vo = 10 eV have been obtained by fitting the energy 
scans and successively they have been kept fixed in the 
recursive fitting procedure. 

The fit quality has been evaluated by means of the 

commonly used reliability factor R — ^ — t~~t ■ 

order to compare our results with previously published 
structural studies, the random error has been estimated 
following the approach of Woodruff and-Bradshaw for 
the analysis of the variable energy PEDnJ Bjj applying 
the method described by Pendry for LEEDJ23 the vari- 
ance of the R-factor is calculated for its minimum value. 
The variance is evaluated as Var{Rjnin) = Rmin\^'^/N, 
where N is the number of independent pieces of strufe 
tural information within the spectrum energy range.tJ 
All the parameter values pi yielding structures with as- 
sociated R-factor lower than Rmin + Var{Rmin) are re- 
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garded as falling within one standard deviation from that 
of the best fit structure. Assuming a parabolic form of 
the R-factor, the random error on the pi parameter is 
evaluated as 



1/2 



2 • Var{R„ 

~5p^ 



(1) 



where the other pj parameters {j ^ i) are taken at 
the minimum of the R-factor. As a result an error of 
±0.03 and ±0.05 A has been found for I and d, respec- 
tively. Due to the intrinsic limit of this approach, that 
assumes the structural parameters to be independent, we 
have also checked, by visual inspection, that a deviation 
of 0.05 A from the best fit value produces significative 
differences in the simulation curve. 

Due to computational limits, the BE and NE energy 
scans have been independently fitted. However, they are 
essentially sensitive to a single parameter (/ and d, re- 
spectively). The preliminary results for I and d, obtained 
by comparing data and simulations by visual inspection, 
have then been used to fit all the polar scans at the same 
time in order to obtain the bond angle a. The new value 
for a has then been used to obtain better I and d val- 
ues from the R-factor analysis of the energy scans. This 
fitting procedure has been iterated until a self-consistent 
evaluation oi I, d and a was obtained for both A and B 
components from the entire set of PED data. 

The refined bond length I and its error have been found 
from the BE data by calculating the reliability R-factor 
for grid simulations on two parameters, ai and c [P = 
af + c^, see Fig. ^) with the other parameters fixed at 
their best fit value. As can be seen in Figs. ^ and|^, a well 
defined minimum is obtained for Ia = 2.82 ± 0.03 A and 
Ib — 2.79 ±0.03 A, while an uncertainty of almost 10° is 
given for the bond direction. We have also observed that, 
in BE conditions, the simulation is weakly affected by the 
position of the surrounding next nearest neighbour Ge 
atoms (not even by the Ge„„„ underneath the Sn atom) . 

In the case of the NE scan, the grid scan simulations 
have been performed spanning over the Sn-Ge„„„ dis- 
tance d and the distance di-n between the first and 
second Ge bi-layers. As shown in Figs. || and ^, the R- 
factor analysis yields two different heights, dA — 2.92 ± 
0.05 A and ds = 3.23 ± 0.05 A, for the A (majority) 
and B (minority) type Sn atoms. It must be emphasized 
that the cluster used for the simulations contains one Sn 
atom only. The relative height between the A and B clus- 
ters is thus undetermined, but the difference — (iliP^ 
- 0.3 A is in, good agreement with the calculatedliO'tSi 
and measurecu'El vertical ripple of the adatom layer. For 
what concerns other structural parameters, the NE scans 
have been found weakly sensitive to d/-//. From a com- 
parison with the polar scan simulations, d/-// has been 
estimated to be 0.05-0.1 A shorter than the Ge bulk 
interlayer separation, for both A and B type Sn atoms. 



The polar scans are sensitive to the whole atomic envi- 
ronment of the Sn atom (particularly to the bond angle 
a), but they display much less features with respect to 
the energy scan (see Fig. The angular patterns have 
thus been fitted to check the consistency of the overall 
structure, as obtained by the NE and BE scans. The 
same bond angle a = 50° has been determined for both 
A and B type Sn atoms, with an uncertainty of ±3°. 

The parameters obtained from the analysis of the 
whole PED datapSet-jare shown in Table | and xampared 
to the calculatedtZrEa and previously measuredOO values. 
There is a qualitative agreement concerning the over- 
all picture. In particular, the nearest neighbour Ge„„ 
atoms nearly follow the Sn atom in its vertical displace- 
ment. According to our data, the three nearest neighbour 
Ge„„ atoms have the same bond length and direction for 
both A and B type Sn atom, within ±0.03 A and ±3° 
respectively. This is also consistent with the ripple of 
the Ge atoms of the firsl||-|lwe£. as reported in the lit- 
erature (- IL|L.- 0.2 A).EH2lO'y The most recent DFT 
calculationsEjiij are in agreement with our results, par- 
ticularly for what concerns the vertical height d. Small 
discrepancies are only obtained for the Sn-Ge„„ bond 
angle, which, on the other hand, yields the higher un- 
certainty in our measurements. Therefore the PED of 
the two components of tjhe-Sn -id spectrum, characteris- 
tic of the (3 X 3) phase^'til yields a|-vertically distorted 
structure, consistently with the XRDeI and XRD/LEEDQ 
determinations. The component with the highest bind- 
ing energy (minority B) is associated to the Sn atoms 
that protrude above the surface. 



IV. CONCLUSIONS 

By comparing the RHEED patterns and the photoe- 
mission spectra we have checked that the correct line- 
shape of the Sn M spectrum on the (3 x,3l phase is 
similar to that reported in several articles. ETEI In addi- 
tion, the lineshape., recently proposed by Kidd et al. for 
the (3 X 3) phasejlj has been only observed when residual 
(2 X 2) phase was present. We have determined the struc- 
ture of the low temperature (3 x 3)-Sn/Ge(lll) phase by 
means of photoelectron diffraction from the Sn 4c? core 
levels. The PED measurements and analysis have been 
performed on a set of photoemission spectra, which can 
be always fitted to two components. Exploiting the scat- 
tering geometry of the experimental apparatus, we have 
separately measured the Sn-Ge„„ bond parameters and 
the vertical height of Sn above the underneath Ge„„„ 
atom for the two kind of Sn adatoms. The results are 
consistent with the mostjECcent XRD measurementsO and 
theoretical predictionsEJ'ES, where one Sn atom out of 
three is vertically displaced by ~ 0.3 A. This vertical rip- 
ple is strongly coupled to the distortion of the first Ge 
bi-layer. In particular, the tetrahedron formed by the 
Sn atom and its three nearest neighbours Ge„„ atoms 
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is found to retain the same structure for both A and B 
components within the experimental error. The vertical 
height above the Ge„„n atom is the main difference be- 
tween the two bonding geometries. The Sn atoms that 
protrude out of the surface have the highest binding en- 
ergy (minority component), while the majority A com- 
ponent (lowest binding energy) is associated to the Sn 
atoms that are pushed closer to the Ge(lll) surface. 
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TABLE I. The best fit values for the parameters a, I, d and di-n, compared to recent experimental measurements and 
theoretical calculations. The calculated random errors are Al = 0.03 A, Ad = 0.05 A, Adi-n = 0.12 A, and Aq = 3° for both 
A and B components. The vertical ripple ASn has been simply evaluated as ds — cJa. The geometric parameters reported for 
the majority A component from Ref. [5] represent the mean value of the two lower height Sn atoms (whose coordinates were 
separately determined) . 
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(V3xV3)30-Sn/Ge( 111) (3x3)Sn/Ge( 111) 



FIG. 1. RHEED patterns (15 kcV) taken from the clean 
surface at Room Temperature (panel o) along the < 112 > 
direction, spots arise from the c(2 x 8) structure; in panel b 
the pattern is taken at Ts = 500 K after about 4/5 of the 
total evaporation time, the (\/3 x \/3) spots are clearly vis- 
ible, while the half-integer ones are very weak and streaky; 
the evaporation is stopped in correspondence of the panel c 
(Ts = 500 K), i.e. when the half- integer peaks completly 
disappear; the RHEED pattern taken along the < 110 > di- 
rection displays the (3 x 3) symmetry upon cooling down to 
150 K (panel d). 
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FIG. 2. The Sn 4d photoemission spectrum taken at 
a photon energy of ~ 210 eV and an overall resolution of 
220 meV, normal emission, TM polarization. The experimen- 
tal data (open circles o) are shown with the fit (full line) to 
four gaussians and a quadratic background (dotted line) . The 
shaded curves indicate the two A and B components whose 
area has been used to determine each data point of the pho- 
toelectron diffraction patterns. Also shown is the spectrum 
taken at a lower coverage (filled circles •), close to the de- 
position stage of panel b in Fig. 1), where residual (2 x 2) 
domains are still present. This spectrum has been taken with 
an overall energy resolution of 120 meV. 




FIG. 3. Upper panel: a sketch of the (3 x 3)-Sn/Ge(lll) 
structure is shown, top view. The underlying matrix of 
small filled circles represents the lattice of nearest neigh- 
bour Genn atoms. The Sn atoms are in T4 absorption 
sites. The unit cells corresponding to the (1 x l)-Ge(lll), 
iVS X V3)R30°-Sn/Ge(lll) and (3 x 3)-Sn/Ge(lll) lattices 
are indicated by the shaded areas. The main symmetry di- 
rections are also indicated. Lower panel: a sideview of the Sn 
atom and its closest Ge atoms. The main fitting parameters 
a, I, d and di-n are also indicated. 
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FIG. 4. Photoelectron diffraction polar scans taken along 
the p"2] direct ion at four different photon energies: 220, 244, 
260 and 360 eV. The emission angle is measured from the 
surface normal (0°). An offset of 0.2 has been added to the 
A x-function for the sake of clarity. At 50-52°, the scans 
display a maximum (220 and 360 eV) or a minimum (244 
and 260 eV) depending on the interference conditions of the 
primary photoelectron wave with that backscattered by the 
nearest neighbour Ge atoms. The best fit simulations are 
also shown (full lines). The vertical dashed lines indicate the 
Sn-Ge„„ bond angle a from the surface normal. 



FIG. 5. Upper panel: X"function of the energy scan in 
bond emission geometry for the A and B components (filled 
and open circles, respectively). This FED measurement is 
mainly sensitive to the Sn-Ge„,j bond length I. Lower panel: 
X-function of the energy scan in normal emission geometry 
for the A and B components (filled and open circles, respec- 
tively). This FED measurement is mainly sensitive to the 
Sn vertical height d above the Gennn atom underneath. Note 
that the data presented here are based on a mcpe extended set 
of spectra with respect to our previous study,L3 thus yielding 
a better statistics and, in the case of normal emission data, a 
more extended energy range too. 
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FIG. 6. Bond length analysis of the main component in bond emission (BE) geometry. Both the photon polarization 
vector and the electron analyzer have been oriented at 50° from the surface normal, along the Sn-Genn bond direction. Left 
upper panel: kinetic momentum (wavenumber) dependence of the x-function for the A (majority) component of the Sn 4d 
photoelectrons in bond emission. The best fit (full line) is overimposed to the experimental data (filled circles •). Right panel: 
contour plot of the MSCD simulation R- factor for the ai and c parameters (see Fig. 3). The bond length I a (dashed thick lines 
with markers on white background) and the bond angle qa (full thick lines with markers on gray background) arc reported for 
a few values. Left bottom panel: statistical analysis of the random error for the parameter Ia with the other parameters at 
their best fit value, as reported in Table 1. The R- factor values (filled squares) have been fitted to a parabolic curve (full line). 
The R-factor variance and the corresponding error on the parameter are also reported. 
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FIG. 8. Vertical height analysis of the main component in normal emission (NE) geometry. Measurements taken in 
transverse-magnetic polarization and with the analyzer oriented along the surface normal. Left upper panel: kinetic momentum 
(wavenumber) dependence of the x-function for the A (majority) component of the Sn 4d photoelectrons in normal emission. 
The best fit (full line) is overimposed to the experimental data (filled circles •). Right panel: contour plot of the MSCD 
simulation R-factor for cIa and dj-u parameters (see Fig. 3); the R-factor values are reported on a white foreground. Left 
bottom panel: statistical analysis of the random error for the parameter dA with the other parameters at their best fit value, as 
reported in Table 1. The R-factor values (filled squares) have been fitted to a parabolic curve (full line). The R-factor variance 
and the corresponding error on the parameter are also reported. 




FIG. 9. Vertical height analysis of the secondary component in normal emission (NE) geometry. The labeling is the same 
of Fig. 8. 
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